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The  purpose  of  this  contract  was  to  expand  the  epicyclic  gear  dynamics 
program  to  add  the  option  of  evaluating  the  tooth,  pair  dynamics  for  two 
epicyclic  gear  stages  with  peripheral  components.  The  option  was 
developed  for  either  stage  to  be  a  basic  planetary,  star,  single  exter¬ 
nal-external  mesh,  or  single  external -internal  mesh.  The  two  stage 
system  includes  an  input  mass  and  shaft,  an  output  mass  and  shaft,  and 
a  connecting  shaft,  where  the  shafts  are  each  modeled  with  torsional 
springs  and  dampers.  The  solution  procedure  was  nearly  the  same  as  the 
procedure  previously  use!  for  determining  tooth  pair  displacements  and 
stresses  in  single  stages.  The  primary  differences  were  that  the  indi¬ 
vidual  gear  component  displacements  were  calculated  rather  than  the  net 
sun-planet  or  ring-planet  displacements.  This  was  necessary  in  order 
to  determine  the  relative  displacements  between  the  shafts  and  the 
input  and  output  gears.  This  generally  increases  the  number  of  degrees 
of  freedom  to  be  solved  per  stage  compared  to  the  single  stage  solu¬ 
tion. 

The  option  to  evaluate  two  stages  makes  the  user's  job  more  difficult. 
Two  stages  of  basic  gear  system  information  must  be  input  as  well  as 
the  additional  component  information.  In  addition,  the  boundary  cordi- 
tions  and  associated  iteration  procedure  become  more  complex.  This  is 
due  to  both  the  increased  number  of  components  and  to  the  time  for  a 
complete  mesh  generally  being  different  for  each  stage. 

A  brief  investigation  into  methods  of  reducing  the  program's  computa¬ 
tion  time  was  done.  The  efforts  focused  on  reducing  the  number  of 
iterations  required  for  boundary  condition  convergence.  It  was  recom¬ 
mended  that  the  beginning  and  ending  values  be  weighted  differently  to 
utilize  previous  iterations  more  effectively.  As  over  90%  of  the  com¬ 
putation  time  is  in  the  numerical  integration  routines,  reducing  the 
number  of  solution  integrations  would  yield  a  direct,  linear  reduction 
in  time. 

Execution  of  the  initial  test  case  indicated  an  instability  in  the 
solution.  The  tooth  pair  load  pattern  is  reasonable-  however  the  mag¬ 
nitudes  of  the  tooth  pair  loads  grow  to  excessive,  \  _  ea? lstic  values 
as  a  function  of  time.  This  could  be  due  to  die  ini.  conditions,  a 
code  error,  or  some  type  of  numerical  instability.  A  procedure  has 
been  recommended  for  eliminating  possibilities  and  determining  where 
the  problem  lies. 


Over  the  past  several  years  -  NASA  nas  developed  a  gear  dynamic  analysis 
and  computer  code  for  standard  and  high  contact  ratio  gears.  The  ana¬ 
lysis  was  expanded  to  include  internal  involute  tooth  forms  in  addition 
to  external  involute  tooth  forms,  then  expanded  to  include  several  pla¬ 
nets  in  epicyolic  gear  systems.  The  development  continued  with  the 
addition  of  helioal  and  double  helical  gears,  a  floating  sun  gear 
option,  a  natural  frequency  evaluation,  a  refined  helioal  gear  com¬ 
pliance  routine,  and  a  flexible  carrier  evaluation. 

The  initial  program  was  developed  for  a  single  spur  gear  mesh  and  to 
operate  over  a  wide  range  of  contact  ratios  (up  to  4.0)  for  analysis  of 
high  contact  ratio  as  well  as  low  contact  ratio  gearing.  This  single 
mesh  program  was  an  extension,  by  Cornell  and  Westervelt  (references  1 
and  2),  of  the  basic  concept  developed  by  Richardson  in  1958  (reference 
3).  The  tooth  pair  compliance  and  stress  sensitivity  formulation  of  the 
single  spur  gear  mesh  program  was  used  in  the  epicyolic  gear  dynamics 
program,  applied  to  each  mesh. 

This  contract,  NAS3-25281,  developed  an  option  for  two  spur  gear  stages 
with  input  and  output  shafts  with  attached  masses  and  a  shaft  connect¬ 
ing  the  two  stages.  This  was  a  practical  extension  to  the  program  as 
more  than  one  gearing  stage  are  often  used  for  speed  reduction,  space, 
weight,  and/or  auxiliary  units.  Thus,  this  extension  allows  for  model¬ 
ing  of  peripheral  components  in  order  to  assess  the  impact  on  the  gear 
tooth  stressing.  Although  the  basic  code  was  completed,  the  cause  of 
the  unexpected  instability  in  the  results  could  not  be  determined  with 
the  cost  constraints  of  the  contract.  NASA,  therefore,  decided  to 
accept  delivery  of  the  code  and  documentation  rather  than  extend  the 
contract  at  this  time. 


TIT. 


The  first  two  tasks  of  contract  NAS3-25821  were  to  modify  the  existing 
Multi -Mesh  Gear  Dynamic  Analysis  Code  (GEARDTNMULT  in  NASA  notation, 
F178  in  Hamilton  Standard  notation)  to  include  a  second  stage  and  per¬ 
ipheral  components .  These  tasks  Included  the  technical  development  and 
roost  of  the  coding  efforts.  Another  smaller  task  was  to  investigate 
methods  of  reducing  computation  time. 


A.  Analytical  Model 


The  dynamic  response  equations  for  a  single  gear  stage  were  expanded  to 
include  additional  degrees  of  freedom  for  an  input  shaft  and  attached 
input  mass,  and  an  output  shaft  with  an  output  mass.  This  was  accom¬ 
plished  by  adding  two  new  equations  and  modifying  other  equations.  Fig¬ 
ure  1  shows  a  diagram  of  the  two  stage  system  with  no  details  for  the 
individual  gear  stages.  Figure  2  shows  a  general  schematic  for  either 
epicyclic  gear  stage.  The  specifics  of  the  gear  tooth  dynamics  and 
program  capabilities  for  each  stage  are  discussed  in  References  1-7. 

For  the  numerical  model,  the  input  and  uu-cput  shafts  each  required  an 
additional  second  order  differential  equation  to  account  for  the  iner¬ 
tias  on  the  input  and  output  shafts. 

The  equation  for  the  input  shaft  and  attached  mass  is: 


(1) 


where  gear  for  stage  1 

hase  radius  of  the  input  gear  to  stage  1 


The  equation  for  the  output  shaft  and  attached  mass  is: 


J  ,  0  ,  ,  +  o  ,  (0  ,  —  0  , ..) 

out  outsh  out  outsh  out2 

+  kout  ^outsh  0out2^  -  Tout  ^ 

The  addition  of  these  equations  changes  the  way  the  torque  is  fed  into 
and  out  of  the  individual  stages.  Previously,  there  was  a  constant 
input  torque  term  in  the  s>on  gear  equation  and  a  constant  output  torque 
term  in  either  the  ring  or  carrier  equation.  The  torque  is  now  fed 
into  the  input  gear  equations  by  way  of  shaft  rotational  displacements 
and  velocities,  which  vary  with  time.  The  rotational  displacements  are 
transformed  into  equivalent  displacements  along  the  line  of  action  for 
each  mesh  of  the  input  gear.  The  output  gear  equations  also  have  the 
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equivalent  shaft  displacement  and  velocity  terms  in  plaoe  of  the  output 
torque  terms. 


The  new  two  stage  program  has  the  capability  of  solving  for  the  dynamic 
tooth  loads  in  planetary  systems,  stair  systems  ,  single  mesh  external - 
external  involute  tooth  forms,  and/or  single  mesh  external -internal 
involute  systems.  Table  1  summarizes  the  available  spur  gear  systems 
for  two  stage  dynamic  solutions.  The  two  stage  program  has  been 
expanded  to  allow  the  user  to  specify  the  input  and  output  gears  The 
single  stage  solution  assumed  the  systems  were  speed  reduction  systems, 
but  the  two  stage  solution  also  allows  for  speed  increasing  systems 
This  option  will  allow  the  user  to  model  a  back-to-back  test  rig,  which 
is  generally  one  reduction  system  and  a  second  similar  system  which 
increases  the  speed.  For  example,  a  single  stage  planetary  system  nor¬ 
mally  has  the  sun  gear  as  the  input  and  the  planet  carrier  for  the  out¬ 
put.  For  the  two  stage  solution,  the  planetary  system  can  have  the 
’-put  component  as  either  the  sun  gear  or  the  carrier  with  the  output 
being  either  the  carrier  or  the  sun  gear  respectively. 

The  equations  for  each  separate  gear  component  were  the  same  as  pre¬ 
viously  developed.  The  basic  equation  development  was  documented  in 
References  4-7,  and  the  equations  are  summarized  in  Appendix  A  for 
convenience.  The  input  and  output  gears  for  each  stage  will  depend  on 
the  specified  system  type,  see  Table  1.  Thus,  the  additional  terms 
needed  in  the  gear  equations,  due  to  the  shafting,  are  written  in  terms 
of  general  input  and  output  gear  equations,  and  are  as  follows. 


input  gear  equation  for  stage  1 

*  cin  (tnl  '  +  kin  (Ytal  '  Yinsh>  "  0 

output  gear  equation  for  stage  1 

*  C12  '  W  +  k12  (W  -  W  '  0 


(3) 

(4) 


input  gear  equation  for  stage  2 
C12  (Yln2  “  Youtl^  +  k12  (yin2  ~  YoutP  "  0 


(5) 


Cout 


output  gear  equation  for  stage  2 
^Yout2  Ycrutsh^  +  ^out  ^Yout2  Youtsh^  ~  ® 


(6) 


where : 


^insh  input  gear  to  stage  *  ^  ^insh 

Youtsh  -  (Rb.  output  gear  from  stage  2^  6outsh 
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Transformations  from  shaft  rotations  to  displacements  along  the  appro¬ 
priate  line  of  action  are  necessary  to  determine  the  tooth  pair  dis¬ 
placements  along  the  line  of  action.  This  allows  the  tooth  pair  dis¬ 
placements  ard  loads  to  be  determined  by  a  procedure  similar  to  the 
single  stage  procedure. 

The  numerical  solution  of  the  dynamic  equations  utilizes  IMSL  (Interna¬ 
tional  Mathematic  and  Scientific  Library)  numerical  integration  rou¬ 
tines.  These  routines  solve  a  system  of  first  order  differential  equa¬ 
tions.  Appendix  B  shows  the  reduction  of  the  second  order  equations  to 
first  order  equations. 

The  tooth  pair  displacements  along  the  line  of  action  are  calculated 
from  the  gear  component  displacements  via: 

v  =  v  -  y  -  y  .  (7) 

isp  2S  ■‘c  Ipi 

v  =  v  .  -  v  -  y  (8) 

*rp  -•pi  ic 

For  the  single  stage  solution,  the  generalized  coordinates  to  be  solved 
by  the  numerical  solution  were  the  tooth  pair  displacements  along  the 
line  of  action  (y  and  y^),  which  minimized  the  number  of  equations 

to  be  numerically  solved,  tor  the  two  stage  system  of  equations,  the 
individual  gear  component  equations  must  be  solved.  This  was  necessary 
in  order  to  obtain  the  input  and  output  gear  displacements  (v  ,  y  ,  y  , 

y  j )  for  both  stages .  Therefore  the  relative  displacement  terms  bet¬ 
ween  the  shafts  and  the  input  or  output  gears  could  be  calculated. 

This  means  there  are  more  degrees  of  freedom  to  be  solved  for  each 
stage  in  the  two  stage  solution  .  There  are  (3  +  number  of  planets) 
degrees  of  freedom  for  each  stage  plus  two  degrees  of  freedom  for  the 
input  and  output  shafts.  Thus,  there  are  a  total  of  (8  +  +  N^) 

second  order  eq1  rations  and  2(8  +  +  N_)  first  order  equations  to  be 

solved  simultaneously. 

The  input  ard  output  shaft  damping  terms  are  calculated  from  a  damping 
ratio,  stiffness,  and  inertias.  The  damping  coefficients  are  then  cal¬ 
culated  using  the  following  formulations. 


2  ?  ll-/  (  1/Jln  *  1  • 

2  5  [>W  (  1,Jout  +  I" 

2  5  C,/  (  W„„t,  * 


B .  Solution  Procedure 


The  solution  procedure  used  to  determine  the  two  stage  dynamic  gear 
tooth  loads  and  stresses  was  very  similar  to  the  solution  for  the 
single  stage  dynamic  gear  tooth  loads  and  stresses.  To  minimize  the 
program  modif ications ,  as  many  of  the  existing  methods  and  procedures 
as  possible  were  used.  Host  sections  of  the  program  needed  some  modi¬ 
fication  in  order  to  process  the  additional  stage.  Many  of  the 
slightly  modified  routines  were  used  twice,  once  for  each  stage,  by 
sending  the  neoessary  information  for  the  relevant  stage. 

The  single  stage  solution  procedure  will  be  reviewed  first,  because  the 
two  stage  solution  procedure  is  essentially  an  extension  of  the  single 
stage  procedure.  The  general  program  procedure  is  illustrated  via  a 
flowchart  in  Figure  3.  In  brief,  the  program  calculates  tooth  pair 
loads  and  solves  the  differential  equations  for  each  of  100  time 
steps — a  piecewise  linear  solution.  The  boundary  conditions  are  deter¬ 
mined  via  an  iterative  procedure,  which  will  be  explained  in  more 
detail  later. 

The  boundary  condition  iteration  procedure  is  intended  to  lead  to  a 
steady  state  condition.  For  the  single  mesh  this  is  accomplished  by 
evaluating  one  full  mesh  cycle,  where  a  mesh  cycle  means  the  time 
required  for  a  tooth  pair  to  move  the  length  of  the  line  of  action. 

The  mesh  cycle  time  for  spur  gears  is  divided  into  100  time  steps  for 
the  numerical  time  stepping  solution,  where  each  increment  of  time  is 
treated  as  linear  with  respect  to  the  tooth  pair  stiffnesses.  The 
boundary  conditions,  for  sun-planet  and  ring-planet  tooth  pair  dis¬ 
placements  and  velocities,  are  compared  at  the  beginning  and  end  of  the 
mesh  time  cycle.  If  the  displacements  and  velocities  are  within  a  spe¬ 
cified  tolerance,  the  system  has  reached  a  steady  state.  If  the  condi¬ 
tions  are  not  within  the  ■colcrcnoe  the  initial  and  final  values  of  the 
displacements  and  velocities  are  averaged  to  obtain  a  new  set  of  ini¬ 
tial  conditions,  and  the  procedure  is  repeated  with  new  initial  condi¬ 
tions.  This  procedure  is  repeated  until  the  boundary  conditions  are 
all  within  the  specified  tolerance.  The  final  iterated  100  time  step 
solution  for  tooth  pair  leads  and  displacements  is  then  post-processed 
for  tooth  stressing. 

The  boundary  conditions  for  two  stages  become  quite  complex  compared  to 
one  stage.  The  main  complexity  is  due  to  the  two  stages  having  differ¬ 
ent  mesh  time  cycles.  This  means  that  steady  state  is  no  longer 
achieved  for  one  mesh  cycle,  rather  it  must  be  a  mutually  common  time 
for  both  stages,  analogous  to  a  common  denominator.  Figure  4  schemati¬ 
cally  shows  two  examples.  The  first  one  shows  a  case  where  the  two 
stages  have  time  cycles  that  can  readily  achieve  a  steady  state,  l.e. 
one  stage  has  a  time  cycle  exactly  2  times  the  other  stage's  time 
cycle.  The  other  example  shows  a  oase  where  47  cycles  would  have  to  be 
evaluated  for  steady  state,  which  would  consume  large  quantities  of 
computer  time  (and  a  code  change  to  increase  the  size  of  an  array). 
However,  it  is  thought  that  this  oase  could  be  approximated  by  evaluat¬ 
ing  five  cycles.  It  follows  that  there  are  systems  in  which  there  is 
no  true  steady  state. 
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Once  the  bov.  ,■  conditions  have  converged,  the  tooth  pair  steady 
state  dlsplcujr .  *2nts  are  post -processed  f  or  stressing  through  the  100 
tine  steps.  If  any  tooth  pair  errors  have  been  requested  for  evalua¬ 
tion.  the  program  continues  to  evaluate  10  tire  cycles  to  simulate  the 
:ser  input  tooth  pair  errors  going'  through  a  mesh  for  either  one  or  two 
stages .  The  method  used  to  process  the  tooth  pair  errors  starts  with 
the  converged  boundary  conditions  from  the  no  tooth  errors  solution  as 
the  initial  conditions  for  the  numerical  solution  with  errors.  The  pro 
gram  thus  simulates  a  transient  response  < for  a  duration  of  ten  mesh 
cycles'  as  the  various  tooth  pair  errors  ccu;e  into  mesh.  The  dynamic 
effects  of  the  tooth  parr  errors  have  generally  been  observed,  tc  dampen 
rut  after  evaluation  of  3-5  mesh  cycles. 

?.::r  the  two  stage  expansion,  it  was  desirable  to  work  with  the  current 
procedure  to  maintain  ccr  morality .  This  implied  that  the  100  time 
steps,  corresponding  to  a  mesh  cycle  tire,  needed  to  be  maintained,  as 
this  fundamental  approach  is  use!  throughout  the  program.  The  first 
question  was  which  stage's  time  cycle  should  be  used  to  'Correspond  to 
tie  ICC  time  steps  embedded  in  the  program.  A  smaller  time  step  for 
the  numerical  solution  would  tend  to  make  for  a  more  stable  numerical 
solution,  while  the  larger  time  step  would  reduce  the  computation  time. 
The  smaller  time  step  was  chosen  for  stabi  11 ty  reasons . 

The  other  question  was  how  to  handle  the  two  different  mesh  times  in 
order  to  achieve  a  "steady  state"  in  a  practical  amount  of  computation 
time.  It  was  decided  that  the  program  could  be  set  up  to  evaluate  the 
iourdary  conditions  after  a  user  specified  number  of  time  cycles  was 
evaluated.  However,  the  fact  that  there  could  be  cases  with  no  true 
steady  state  made  it  desirable  to  try  and  approximate  a  steady  state 
and  then  evaluate  a  transient  response  using  the  approximated  boundary 
conditions  as  initial  conditions.  It  is  possible  that  the  instability 
could  be  related  to  these  approximate  steady  state  boundary  conditions. 

In  order  to  verify  the  approximate  steady  state  boundary  condition  con  ¬ 
cert,  the  original  single  mesh  program  was  used  to  determine  if  a  small  1 
charge  in  the  mesh  tire  cycle  would  significantly  affect  'he  converged 
bovidary  condition  values.  Several  cases  were  run  with  slightly  dif¬ 
ferent  numbers  of  teeth  (which  determines  the  mesh  time).  The  iterated 
results  showed  minimal  variation  in  the  converged  values  for  displace¬ 
ment  and  velocity.  Thus,  this  implied  that  it  would  be  possible  to 
evaluate  'approximate  two  stage  boundary  conditions'  at  a  time  when  the 
two  stages  were  close  to  a  steady  state.  The  program  will  use  the  ini¬ 
tial  conditions  obtained  from  the  approximate  steady  state  solution  and 
continue  to  evaluate  10  additional  cycles  for  a  transient  response. 

The  transient  response  utilizes  the  same  procedure  used  to  evaluate 
tooth  pair  errors  in  the  single  stage  code.  The  solution  procedure, 
with  IX  time  steps  per  cycle,  is  essentially  the  same  except  there  is 
no  boundary  condition  averaging .  It  should  be  noted  that  the  two  stage 
solution  compares  the  tooth  palp  displacements  and  velocities  for  con¬ 
vergence  .  not  the  individual  gear  components.  The  initial  aid  final 
oorxiitions  for  the  individual  components  fare  each  averaged  aid  that 
average  is  used  for  the  initial  conditions  for  the  next  iteration. 


The  subroutines  written  for  the  two  stage  solution  generally  follow 
similar  station  to  that  used  in  previous  codes,  both  in  subroutine 
names  and  variable  names.  The  code  was  written  assuming  a  itexbnuxn  of  10 
planets  pen  stage,  spur  gears  only,  and  no  floating  sun  or  flexible 
carrier  option. 

Several  new  routines  and  some  slightly  modified  routines  were  written 
in  order  to  accommodate  the  additional  information  required  for  two 
different  stages.  The  number  of  planets  per  stage  Is  limited  to 
10.  which  should  be  more  than  sufficient  for  any  real  system.  This  is 
a  reduction  in  the  number  of  planets  allowed  in  the  single  stage 
solution,  but  simplifies  the  code  changes  by  utilizing  the  existing 
arrays  and  the  corresponding  dimensions.  The  arrays  now  contain  stage 
i  information  in  the  first  10  elements  and  stage  2  information  in  the 
second  10  elements ,  for  the  same  total  number  of  array  elements  as  the 
previous  solution. 

Hie  procedures  and  calculations  of  the  new  program  are  essentially  the 
same  as  the  single  stage  program.  The  primary  charges  come  through 
nearly  every  section  or  subroutine  being  executed  for  stage  1  and  then 
for  stage  2.  Some  of  the  basic  parameters,  such  as  mesh  time  cycle, 
required  additional  logic  in  order  to  implement  the  option  of 
reversed"  systems  (where  the  systems  are  speeu  increasers  rather  than 
reducers).  Most  of  the  minor  modifications  were  made  to  accommodate 
the  additional  parameters  of  the  second  stage,  and  did  not  change  any 
basic  concepts.  Because  of  the  option  of  multiple  cycles  being  evalu¬ 
ated  for  boundary  condition  convergence,  it  was  also  necessary  to  add  a 
dummy  time  parameter  to  obtain  proper  tooth  pair  contact.  Without  the 
dummy  time  parameter,  the  stage  with  the  shorter  mesh  time  cycle  lost 
all  tooth  pair  contact  during  the  second  100  time  steps.  Thus,  to 
obtain  proper  contact  evaluation,  a  dummy  time  was  reset  to  zero  at  the 
beginning  of  each  mesh  time  cycle  for  the  stage  with  the  smaller  mesh 
time  cycle. 

The  input  routines  were  modified  to  meet  NASA  programming  standards. 

The  additional  stage  of  information  necessary  for  the  two  stage  solu¬ 
tion  required  modification  of  the  code  to  set  up  the  input  variables. 
This  portion  of  the  code  was  previously  in  the  preprocessing  routine; 
however,  the  preprocessing  routine  (REAB2)  is  a  very  large  routine  and 
does  not  meet  the  150  executable  line  requirement.  Thus,  the  section 
of  the  code  for  setting  input  variables  to  meaningful  notation  was 
modified  to  be  separate  subroutines,  accepting  input  for  either  one  or 
two  stages  of  Information,  and  the  remainder  of  READ2  remained 
■uncharged  except  for  additional  arguments.  Three  subroutines  were 
necessary  to  convert  a  section  of  the  old  routine  into  acceptable 
length  routines  and  a  fourth  was  necessary  to  control  which  stage  was 
being  preprocessed. 
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this  task  was  to  emphasize  the  reduction  in  computer 
time  required  it/  the  program.  In  analyzing  the  present  CPU  usage  of 
the  program,  it  was  determine!  that  o*rer  90  percent  of  the  computa¬ 
tional  time  was  in  tiie  numerical  integration  routine.  The  basic  alter¬ 
natives  to  reducing  program  solution  tuae  were  as  follows. 


Reduce  the  number  of  time  steps  used  throtjgh  the 

time. 
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Item  a.  v.cis  thought  to  be  practical  for  single  mesh  cases  only,  since 
:.:dy  .tit  these  oases  could  the  timing  of  events  through  the  mesh  time 
cycle,  be  accurately  estimated,  item.  b.  was  not  thoroughly  explored, 
but  was  felt  to  hold  some  promise .  Item  c.  was  the  main  focus  of  this 
study,  and  :x.r.:i.i  be  applicable  to  any  level  of  the  code,  i.e. ,  single 
mesh ,  mrtltiple  mesh,  arid  multiple  stages. 

In  crier  to  investigate  better  methods  of  convergence,  the  flexible 
carrier  test  -case  given  in  Reference  4  was  used.  This  case  was  ini¬ 
tially  solved  by  calculating  input  conditions  as  specified  in  Appendix 
A  of  Reference  4,  and  marry  boundary  condition  Iterations,  which  use 

significant  CPU  time. 

Figure  5  is  a  plot  of  the  ending  values  of  the  sun  gear-  displacement 
and  velcoity  versus  the  number  of  iterations  through  a  gear  mesh 
yeri  d .  1 1  shews  that  the  method  of  selecting  new  values  is  creating  a 

cyclic  behavior  that  will  take  many  iterations  to  converge.  The  mean 
value  of  the  plotted  functions  match  the  converged  value  of  the  solu- 
t:  or..  The  behavior  of  the  planet  axd  ring  gear  results  was  similar. 

In  crier  t  ;  evaluate  the  convergence  duf ficul ties .  the  code  was  tempo- 
rar:  ly  luugoi  to  weigh  the  beginning  and  ending  values  differently  in 
r  r.ier  to  obtain  a  higher  weighting  to  previous  iterations.  The  revised 
code ,  shown  in  Figure  6,  has  not  been  incorporated  in  the  new  program, 
hut  was  used  for  this  investigation  only.  The  conv..rgence  could  be 
tailored  for  a  particular  case  by  changing  coefficients  of  the  dls- 
plir*  inerts  from  the  previous  iteration  and  the  current  iteration  (e.g. 
.TTA  and  X30  respectively.; ,  and  similars/  for  the  velocities  (SXS1  and 
>17 1  .  The  sum.  of  tie  ooeff icier, is  must  he  equal  to  unity.  The  results 
of  -using  j  7  and  0.3  for  the  ui.-placenert ,  and  0.5  and  0.5  for  the 
vc  1  /cl tv  are  shewn  in  Figures  V  aiwi  6 . 


In  order  to  reduoe  the  computation  time  for  a  dynamic  solution,  it  is 
recommended  that  this  method  be  added  to  the  code,  with  the  weighting 
coefficients  as  input  parameters  rather  than  fixed  values  in  the  code. 
It  is  further  recommended  that  code  be  written  to  automate  the  plots  as 
shown  in  Figures  7  and  8.  The  resulting  cede  should  substantially 
enhance  the  user's  understanding  and  the  convergence  of  the  boundary 
conditions. 


D.  Discussion  of  Instability  and  Test  Cases 

Task  Ic  was  to  evaluate  the  two  stage  model  and  determine  the  Interac¬ 
tion  between  two  stages  and  the  effects  of  input  and  output  shafting 
and  attached  masses.  This  was  not  fully  accomplished  due  to  the  ini¬ 
tial  test  case  showing  an  instability  in  the  results.  However,  same 
general  comments  will  be  made  regarding  the  apparent  instability  and 
other  observations  from  a  variety  of  test  oases. 

The  first  test  case  that  was  used  approximated  the  OH-58  dropped  tooth 
planetary  design  being  tested  at  NASA  Lewis.  A  dropped  tooth  design 
"drops"  a  tooth  from  each  of  the  planets  so  that  the  number  of  teeth  on 
a  planet  divided  by  the  number  of  planets  is  not  an  integer  number. 

The  purpose  of  the  design  1s  to  reduoe  dynamio  tooth  loading.  The 
actual  planetary  system  could  not  be  modeled  with  the  program,  because 
the  dropped  tooth  design  leads  to  different  pressure  angles  for  the 
sun-planet  and  ring-planet  meshes.  The  program  assumes  equal  pressure 
angles  for  the  sun-planet  and  ring-planet  mesh,  although  an  equivalent 
buttress  tooth  form  might  be  a  feasible  modeling  method.  Figure  9a 
shows  the  input  test  case  for  a  single  stage  and  Figure  9b  shows  the 
corresponding  output  summary.  This  case  had  assumed  the  pressure  angle 
of  the  sun-planet  mesh  for  both  the  sun-planet  and  ring-planet  meshes 
and  used  the  program's  geometry  preprocessor. 

The  next  step  was  to  evaluate  two  stages,  simulating  the  back-to-back 
test  rig  at  NASA.  An  example  Input  data  set  for  two  stages  is  shown  in 
Figure  10.  The  first  stage  was  modeled  as  a  planetary,  with  the  sun 
gear  as  the  input  gear,  and  the  second  stage  was  a  "reversed  planetary" 
with  the  carrier  as  the  Input.  This  case  was  used  for  initial  verifi¬ 
cation  of  the  code  Interaction  and  debugging  purposes.  The  geometry  for 
the  two  stages  was  the  same,  with  the  primary  difference  between  the 
stages  being  the  Input  torque  and  rpm.  The  different  torques  led  to 
slightly  different  tooth  pair  stiffnesses,  as  the  torque  level  influen¬ 
ces  the  Hertzian  component  of  the  stiffness  function,  see  reference  2. 
The  difference  in  rpm  oould  affect  the  dynamio  load  levels,  particu¬ 
larly  if  near  a  resonant  speed. 

Table  3  summarizes  the  test  oases  run  with  the  two  stage  program  and 
the  resulting  component  displacements.  The  initial  test  oase,  oase  1, 
was  designed  to  uncouple  stage  1  from  stage  2.  Thus,  the  connecting 
shaft  stiffness  was  very  low.  The  input  and  output  shaft  stiffnesses 
were  also  low,  because  some  preliminary  test  case  results  indicated 
large  displacements  for  very  stiff  input  and  output  shafts.  The  solu¬ 
tion  seemed  to  be  very  sensitive  to  the  values  of  the  shaft  stiff- 
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nesses,  as  indicated  by  comparison  of  case  4  and  5  (stage  1).  For  a 
decrease  in  input  shaft  stiffness  by  a  factor  of  ten,  the  magnitudes  of 
the  planetary  component  displacements  along  the  line  of  action 
decreased  by  orders  cf  magnitude . 

Some  of  the  dynamic  load  results  for  the  initial  test,  case  have  been 
plotted  in  Figure  11.  These  results  illustrate  a  typioal  tooth  pair 
meshing  pattern  cf  engagement  and  disengagement  of  the  teeth.  The 
meshing  pattern  also  repeats  for  subsequent  cycles,  as  would  be 
expected.  What  is  unusual  about  these  test  results,  is  the  substantial 
increase  in  the  rraghitude  of  the  loads  as  time  progresses.  This 
increasing  magnitude  is  also  illustrated  in  Figure  12,  which  shows  the 
maximum  stress  in  a  cycle  vs.  the  transient  cycle  number.  It  is  evi¬ 
dent  from  this  figure  that  there  is  some  type  of  instability  in  the 
solution.  It  is  interesting  that  the  Hertz  stress  increases  in  a 
nearly  linear  fashion,  while  the  maximum  bending  stress  is  a  higher 
order  function.  This  is  logical  in  that  a  nonlinear  increase  in  loads 
would  correspond  to  a  nonlinear  increase  in  bending  stress.  It  is 
recammended  that  the  system  equations  be  evaluated  for  stability. 

The  cases  summarized  in  Table  3  show  the  effect  of  varying  the  inprut 
and  output  masses  and  shaft  stiffnesses.  Only  one  boundary  condition 
iteration  was  performed,  which  was  sufficient  to  observe  the  change  in 
the  component  displacements  due  to  the  modified  parameters  ) .  The 
affect  of  varying  the  masses  on  the  input  or  output  shaft  was  inversely 
proportional  to  the  input  or  output  shaft  displacements  respectively. 
For  example,  the  input  mass  for  case  2  is  10  times  the  input  mass  for 
case  1  aid  the  input  shaft  rotational  displacement  decreased  by  the 
same  factor  of  10.  The  other  stage  1  component  displacements ,  sun, 
carrier,  and  planets,  also  showed  decreases  of  nearly  the  same  factor. 
This  indicates  that  the  shaft  displacements  are  dominaling  the  individ¬ 
ual  input  and  output  gear  equations,  and  the  tooth  pair  dynamics  axe 
therefore  overshadowed . 

Case  3  is  the  same  as  Case  1  except  the  input  and  output  masses  have 
been  reversed.  The  input  and  output  shaft  displacements  were  corre- 
spcndingly  interchanged .  As  would  be  expected,  the  gear  displacements 
for  the  two  stages  were  interchanged  also,  but  with  some  variation  due 
t.o  the  change  in  torque  and  rpm  of  the  two  stages. 

In  Case  4,  a  very  stiff  shaft  connecting  the  two  stages  was  investi¬ 
gated.  This  case  verified  that  the  displacement  of  the  output  gear 
from  Stage  1  (planet  carrier)  would  equal  the  displacement  of  the  input 
gear  to  Stage  2  (planet  carrier)  for  a  rigid  connection. 

Case  5  was  a  check  case  to  verify  that  the  input  mass  and  shaft  could 
Ijs  isolated  from  Stage  I  by  using  zero  torsional  shaft  stiffness.  By 
using  a  low  connecting  shaft  stiffness,  this  case  also  verified  that 
the  two  stages  could  be  isolated  from  each  other.  By  having  zero  tor¬ 
si  oral  input  shaft  stiffness,  there  is  no  torque  transmitted  to  the 
input  gear  cf  Stage  1  aid  therefore  approximately  zero  displacements  in 
]  The  seoord  stage  displa.deme.nts  were  nearly  the  same  as  Case 
4 .  which  verified  the  stages  were  isolated  from  each  other. 


Case  6  was  nearly  the  same  as  Case  4,  but  the  connecting  shaft  stiff¬ 
ness  was  a  lower  and  more  realistic  value.  All  of  the  component  dis¬ 
placements  were  slightly  larger  for  the  softer  shaft.  There  was  also  a 
snail  difference  in  the  magnitudes  of  the  output  gear  from  Stage  1  and 
the  input  gear  to  Stage  1.  This  is  logical,  because  as  the  connecting 
shalt  stiffness  decreases,  the  relative  displacement  of  the  two  ends 
will  increase. 

Case  7  increased  the  torsional  stiffness  of  the  input  shaft  by  several 
orders  of  magnitude.  The  resulting  displacements  shown  in  Table  4 
indicate  the  displacements  also  increased  by  orders  of  magnitude  for 
the  first  stage  components.  The  second  stage  component  displacements 
also  increased  substantially,  but  not  as  much  as  the  first  stage. 

There  are  several  possibilities  for  the  sauroe  of  the  instability. 
Figure  13  illustrates  a  procedure  for  eliminating  possibilities  and 
tracking  down  the  problem.  The  recommended  procedure  would  start  with 
a  very  simple  "two  stage"  system,  as  shown  schematically  in  Figure  14. 
This  simple  system,  effectively  a  gear  train,  with  constant  tooth  pair 
stiffnesses  could  be  evaluated  vising  a  small  program,  separate  from  the 
gear  program.  The  system  should  have  a  1:1  gear  ratio  for  simplicity 
and  to  ensure  a  steady  state  can  be  achieved  for  boundary  condition 
convergence.  If  the  simple  system  shows  an  unstable  behavior,  it  will 
be  known  that  there  is  a  problem  with  the  basic  problem  formulation. 
That  is  the  least  likely  souroe  of  the  problem,  but  it  is  the  logical 
starting  point.  The  model  complexity  can  then  be  increased  by  adding 
in  the  variable  tooth  pair  stiffness.  If  the  independent  solution 
yields  reasonable  results,  the  same  system  can  be  run  in  the  full  mul¬ 
tiple  mesh  gear  program.  If  the  results  are  inconsistent,  it  would 
indicate  an  error  in  the  two  stage  multiple  mesh  code. 

If  the  results  are  consistent  for  the  1:1  gear  ratio  oases,  the  pro¬ 
blem  is  most  likely  related  to  the  initial  conditions  and/or  the  boun¬ 
dary  condition  iteration  scheme.  This  would  indicate  a  better  method 
for  estimating  the  initial  conditions  and  for  determining  convergence 
should  be  devised.  It  should  be  noted  that  this  is  not  minor  task.  In 
general,  each  stage  plus  an  input  or  output  shaft  will  have  2  times  (4 
+  number  of  planets)  for  boundary  conditions  to  be  determined.  Not 
only  is  it  difficult  to  obtain  initial  estimates  for  the  tooth  dis¬ 
placements  and  velocities  for  each  component  of  the  epicyclic  stage, 
but  it  is  even  mor°  difficult  to  determine  a  general  method  of  combin¬ 
ing  the  two  stages. 

Another  test  case  was  tried  to  determine  if  the  instability  was  system 
dependant.  The  second  case  chosen  was  the  Stoeckicht  2K-H  planetary 
from  Reference  5.  The  two  stage  test  case  was  set  up  as  two  planetary 
reduction  stages  with  identical  geometry  and  basic  gear  parameters. 

The  only  difference  between  these  two  stages  would  be  the  torque  and 
the  rpm.  This  case  exhibited  similar  instabilities ,  with  tooth  loads 
increasing  substantially  with  time. 
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iv .  CtOUqIiIyI  I 


An  option  to  evaluate  two  gear  stages,  with  a  connecting  shaft  and 
input  and  output  shafts  with  attached  masses,  has  been  added  to  the 
epi cyclic  gear  tooth  dynamic  analysis  program.  There  were  extensive 
code  modifications  to  accommodate  two  stages  of:  geometry  preprocess¬ 
ing,  tooth  pair  load  calculations ,  and  dynamic  equations.  The  new  code 
utilizes  the  same  general  methods  and  procedures  as  were  previously 
used  for  the  single  stage  solution  ard  includes  additional  degrees  of 
freedom  for  the  attached  masses.  The  two  stage  option  also  allows  for 
speed  increasing  systems  as  well  as  speed  reduction  systems. 

It  is  recommended  that  code  changes  be  incorporated  for  a  more  robust 
boundary  condition  iteration  scheme.  The  proposed  code  changes  can 
significantly  reduce  the  number  of  iterations  required  for  a  converged 
steady  state  solution,  and  therefore  reduce  the  computation  time  for  a 
dynamic  solution. 

The  program's  plotting  capabilities  should  be  modified  to  use  more  uni¬ 
versal  plotting  routines.  The  plotting  capabilities  should  also  be 
expanded  to  process  the  large  amount  of  output  from  a  two  stage  solu¬ 
tion.  In  addition,  plotting  capabilities  for  the  boundary  conditions 
should  be  added  to  enhance  the  user's  understanding  of  the  boundary 
condition  convergence  process. 

The  initial  test  case  indicated  an  instability  in  the  dynamic  solution, 
where  the  gear  tooth  displacements  continually  increase  with  time,  bey¬ 
ond  any  reasonable  values.  Several  variations  on  shaft  stiffnesses, 
attached  masses  end  shaft  damping  were  evaluated  to  investigate  the 
problem.  The  tooth  pair  loads  with  respect  to  time  indicated  the  tooth 
pair  contact  patterns  were  reasonable ,  in  that  the  basic  pattern  was 
repeated  with  time,  although  the  magnitudes  continued  to  increase. 
Further  evaluation  of  the  problem  could  not  be  made  because  of  fund 
limitations . 
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Appendix  A.  Spicyclic  Gear  System  Equations  for.  .Either.-Stage 


The  following  second  order  equations  of  motion  describe  the 
individual  gear  components  in  an  epicyclic  system.  They  are 
wit ten  in  terms  of  the  displacements  along  the  lines  of  action. 
The  development  of  the  basic  equations  is  document©!  primarily  in 
Reference  1  and  6. 


Sun  gear  equation: 


N  N 

n  y  ’■Id  v  -t-  I  L  =  r .  /R 

s/s  i=1  SPi^i  i_l  sp±  \ 


(A.l) 


Planet  gear  equation: 


mv-d  v  +  d  v  -  L  +  L  =0 

V Pi  ^PiysPi  rPi  rPi  ^i  rPi 


(A. 2) 


ng  gear  equation: 


m  y 
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(A. 3) 


Carrier  equation: 


m  v 
c  ■‘o 


N  ,  N  ,  N  N 

-Id  v  -Id  v  -IL  -IL 

i-1  rpi  rpi  i-1  1-1  rPi  1-1  **i 


(A.  4) 


NOTE:  For  the  two  stage  solution,  the  constant  torque  terms  indi¬ 
cated  above  go  to  zero.  The  "torque"  is  input  to  the  appropriate 
gear  equations  via  terms  for  relative  displacement  with  respect 
to  a  shaft. 


The  tooth  pair  load  terms  are  determined  from  the  following  equa¬ 
tions  . 


=  y  f(y  -  e  -  X2  02  )ti  <j>  ]  (A. 5) 

>1  SPi  SPii  SPii  SPin  sp^sp,  J 
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L  =  V  [  (y  -  e  -  X2  B2  )n  <p 
rp.  .  , 1  *rp.  rp  . .  rp  . .  rp  . .  rp  . .  rp..- 

\Ji  *Ji  Ji  Ji 


(A. 6) 
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Appendix  B :  L'  rst  Order  Differential  Equations 


The  numerical  integration  routines  solve  for  a  system  of  first 
order  equations.  The  transformation  of  the  second  order  differ¬ 
ential  equations  to  first  order  equations  is  shown  below.  The 
subscripted  numbers  on  the  derivatives  correspond  to  elements  in 
the  array  (XI)  in  the  program. 

Input  shaft  equation  reduction: 


Let 


*3 


Then 


*3 


e.  , 

insh 

*3  ~  9insh 

[  t.  -  C.  (x.  -  0.  . )  -  k.  (x.  -  0.  .)]/J. 

1  in  n  4  ini  in  4  ini  J  in 

*3 


(B.l) 

(B.2) 

(B.3) 

(B.4) 


Output  shaft  equation  reduction: 


Let 

«  •• 

x5  ~  outsh 

*6  =  X5  =  9outsh 


(B.  5) 
(B.6) 


Then 

x5  =  [  Tout  ~  cout^6  "  ®out2-)  ~  kout(x6  "  0out2^/Jout  (B‘7) 
x6  =  x5  (B'8) 


The  remaining  equations  will  vary  depending  on  the  specific  sys¬ 
tem  type,  see  table  of  system  types  in  the  main  report  text.  For 
simplicity,  the  following  first  order  equation  derivations  will 
assume  both  stage  1  and  stage  2  are  planetary  systems.  Thus,  the 
generalized  input  and  output  gear  equations  in  equations  3-6 
are  now  specified,  where  the  input  gear  is  the  sun  gear  and  the 
output  gear  is  the  planet  carrier  for  both  stages. 

Sun  gear  equation  for  stage  1: 


Let 


i  ‘^s  9inl  ^bs 


X2  =  X1 


(B  .9) 
(B.10) 
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Then 


N1  N1  . 

X1  '  '  [  if ldSP±  +i^LsPi  +  C^(X2  "  Vt4) 

+  km(x2  -  hsV  I  '  ”s 


X2  -  X1 


(B. 11) 
(B. 12) 


Carrier  equation  for  stage  1 : 
Let 


Then 


*7  = 

*6  =  *7  =  i 


(B. 13) 
(B. 14) 


N1 

*7  -  <  E  * 


N1  ,  N1  N1 

_ y  +  I  d  y  +  EL  +  IL 

i-1  rPi  rPi  i-1  ^i  i=l  **i  i-1  rPi 


C12(X8  X12^  k12(X8  X12^  /  mc 

x8  =  X? 


CB. 15) 
(B. 16) 


Ring  gear  equation  for  stage  1: 

Xg  -  yr  =  0.0  (B. 17) 

X10  =  ^  =  Yr  =  0.0  (B. 18) 

Note:  ring  gear  is  fixed  for  a  planetary  system. 


Planet 

gear  equations  for  stage  1 : 

Let 

x.  -a  =  y  .  for  i  =  1  to  N1 

i+16  Ipi 

(B. 19) 

x.  =  x,  .e  for  i  =  1  to  N1 

i+Nl+16  i+16 

(B.20) 

Then 

x.  =  f  d  v  -d  v  +L  -L  ]  /  m 

i+16  L  spi*spi  rpixrpi  spi  rpi  Pi 

(B.21) 

• 

Xi+Nl+16  =  Xi+16 

(B.22) 

17 


Sun  gear  equation  for  stage  2: 
Let 


X11 

■  ?s2 

(B.23) 

X12 

‘  X11  •  ?S2 

(B . 24) 

Then 

X11 

N2  .  N2  . 

[i=ldsp2i  *sp2±  +i=1Lsp2i  +  C12(X12  '  V 

+  k12(x12  -  V  1  7  ms2 

(B.25) 

X12 

"  xn 

(B.26) 

Carrier  equation  for  stage  2: 

Let 

X13 

• »  •• 

^c2  ~  0out2  ^102 

(B. 27) 

X14 

=  X13  =  Yc2 

(B.28) 

Then 


,  N2  N2  N2  N2 

X13  '  I^VA-pS.  +  hV  +  JjVj. 

”  GoutrRbc2^6  “  X14^  “  1’:out^Rbc2X6  ~  X14'>  ^  1  mc2 


(B.29) 


x 


14 


=  X 


13 


(B.30) 


Ring  gear  equation  for  stage  2: 

x15  =  0.0  (B.31) 

x16  "  x15  ’  ?r2  ■  0  0  (B'32) 

Note:  ring  gear  is  fixed  for  a  planetary  system. 
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Planet  gear  equations  for  stage  2- 


Let 


v  =  x  .  for  i  =  1  to  N2 

Xi+2N1+16  pi 

Xi+2N1+N2+16  =  Xi+2N1+16 
xi'2Nl-l6  =  ^  dsp2i^sp2i“  drp2i^rp2i+  Lsp2i 

-  Lrp2j  7  Vp 


Xi+2N1+N2+16 


Xi+2N1+16 


(B. 33) 
(B  .34) 


(B . 35) 
(B. 36) 
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Appendix  C. 


C.^  =  damping  coefficient  for  the  input  shaft  (in-lbs-s/rad) 
C„ =  damping  coefficient  for  the  output  shaft  (in-lbs-s/rad.) 
C^0  =  damping  coefficient  for  the  output  shaft  (in-lbs-s/rad) 


=  tooth  pair  damping  (lb-s/in) 

2  =  tooth  spacing  error  (in) 

o  w  .  .  ^ 

*•11 

j.  =  rctatioral  inertia  attached  to  input  shaft  (in-lbs-s  2 /rad) 

=  rotational  inertia  attached  to  output  shaft  (in-lbs-s 2 /rad) 

k.  =  torsional  stiffness  of  the  input  shaft  (in-lb/rad) 
k  ^  -  torsional  stiffness  of  the  output  shaft  (in-lb/rad) 

k.  0  =  torsional  stiffness  of  the  shaft  connecting  the  tvo  stages 

L  *o 

I  =  tooth  pair  leads  for  planet  mesh  i  (lbs) 

-  =  rotational  (equivalent)  mass  (lb-s2/in) 

h  =  number  of  planets 

R  =  base  radius  (in) 

7  =  displacement  along  the  line  of  action  (in) 

£  =•  cam  modification,  ring-planet  mesh 

rn  ■  . 

*  j- 

2  =  cam  modification,  sun-planet  mesh 

*  J-*- 

r  =  ring -planet  tooth  pair  spring  rates 


SU. 
-  1 


'  1 1 T.  Di.SL.  u* 


or  spring  rates 


9  =  rotational  displacement  (rad) 
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’  darning  ratio 

-  input  torque  (in- lb) 

-  output  torque  (in- lb) 

-  identity  function  for  ring-planet  tooth  pair 

-  identity  function  for  sun-planet  tooth  pair 

-  o  or  1  dererding'  on  whether  the  tooth  ccntac 
the  rrofile  r.cdif uonticn  can  or  net 


-•h-'A  1 

-*-t  -  * 

tape  i 

^  i 


T’X-I 


Table  1 :  Spur  Gear  Systems  Available  for  Two  Stage  Dynamic 

Solution 


System 

:unber 

Spur  Gear 

System  Type 

Input  Gear 

Output  Gear 

planetary  system 

sun  gear 

planet  carrier 

planetary  system 

planet  carrier 

sun  gear 

O 

star  system 

sun  gear 

ring  gear 

o 

star  system 

ring  gear 

sun  gear 

4 

single  external- 
external*  mesh 

sun  gear 

planet  gear 

4 

single  external- 
external*  mesh 

planet  gear 

sun  gear 

r 

single  external- 
internal**  mesh 

planet  gear 

ring  gear 

5 

single  external- 
internal**  mesh 

ring  gear 

planet  gear 

*  External-external  mesh 

means  a  pinion  and  a 

gear,  both  with 

external  involute  tooth  forms . 


**  External -internal  mesh  means  a  pinion  and  a  gear,  with  the 
inion  having  an  external  involute  tooth  form  and  the  gear 
,.aving  an  internal  tooth  form. 
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thcj  fU  pJ  c~t"  rf 


Table  2 : 


System  Type:  Planetary 

Sneed  Range:  Up  to  1620  rpm  (sun  gear) 

Torque  Range:  Up  to  12,450  in-lb  (sun  gear) 

Number  of  Planets:  3 

age  1  input :  Sun  gear 

age  2  input:  Planet  Carrier 

ge  1  output:  Planet  Carrier 

ge  2  output ;  Sun  gear 


heral 


Input  Output  Connecting 


Tr.ertia  (in-lb-s  )  2.5  1.0 

-.iff ness  (in-lb/rad)  l.OxlO5  1.8xlOV  3.9xl09 


lane  Cary  Gear 
haracteristic 

Sun 

Gear 

Planet 

Gears 

Ring 

Gear 

Carrier 

umber  cf  Teeth 

27 

35 

99 

- 

uametral  Pitch* 

6.8571 

8.6571 

9.14286 

- 

ressure  Angle* 

(degrees ) 

24.6 

24.6 

20.19 

- 

tot  Radius  (in) 

1 . 39 

1 . 802 

5 . 554 

- 

ooth  Tup  Radius  (in) 

1 .655 

2.067 

5.35 

- 

ace  Width  (in) 

1 . 375 

1 . 20 

1  .  CO 

- 

r.ortia  (in-lb-s*) 

0 . 026 

0.016 

0.0 

C .  20 

The  variation  in  dianu  ral  pitch  and  pressure  angle  cannot  be 
accurately  modeled,  as  +he  program  assumes  equal  pressure  angles 
at  the  sun -planet  and  ring-planet  meshes. 
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Table  3 : 


Two  Stage  Test  Case  Parameter  Variations 


E  NUMBER 

ATTACHED  INERTIAS 

(ln-lbs-s1)  ! 

SHAFT  STIFFNESSES 
(in-lb/rad.) 

INPUT 

OUTPUT 

INPUT 

OUTPUT 

CONNECTING 

1 

10 

15 

10 

| 

10 

10 

2 

100 

10 

10 

10 

10 

3 

15 

10 

10 

10 

10 

4 

10 

10 

10 

10 

lxlO9 

5 

10 

10 

0 

10 

10 

6 

10 

10 

10 

10 

lxlO5 

r* 

10 

10 

lxlO6 

10 

10 
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Table  4: 


DISPLACEMENTS,  IN. 


CASE 

NUMBER 

STAGE 

INPUT 

SHAFT 

xlCT 

OUTPUT 

SHAFT 

xlO* 

SUN 
GEAR, 
xl  04 

CARRIER 

Xl05 

PLANET 

4 

xio* 

1 

1 

.8017 

.6941 

.7872 

.2067 

2 

- 

.5345 

.4274 

.  503 

.  1292 

1 

.08018 

— 

.06955 

.07985 

.02071 

2 

- 

.8018 

.6408 

.7528 

.1937 

3 

1 

.  5345 

— 

.463 

.526 

.  1379 

2 

- 

.8018 

.6408 

.7534 

.  1937 

4 

1 

.8003 

- 

.6899 

.7629 

.2051 

2 

- 

.8003 

.639 

.7629 

.  1928 

5 

1 

.8004 

- 

. 616E-4 

. 1095E-2 

. 1974E-4 

2 

- 

.8004 

.638 

.7466 

.  1925 

6 

1 

.6017 

— 

.6935 

.777 

.2065 

2 

- 

.8017 

.6414 

.764 

.  1938 

7 

l 

19.889 

- 

2755.6 

9612 

1247.9 

2 

13.364 

59.4 

231 . 8 

26.1 

NOTE  : 

These  displacements 

are  at  the 

end  of  the  first  boundary 

condition  iteration. 
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ini 

_ _ i 

Stage 


outl 

1 


w±n2 


wout2 


Stage  2 


outsh 


Figure  1 :  Two  Stage  System  Diagram 


Read  input  data 
for  1  or  2  stages 
and  store  as  needed 


Preprocess  each  stage 
for  geometry  S'  other 
program  variables 


Call  appropriate  routine 
for  calculating  dynamic 
equation  constants  for 
either  1  or  2  stages 


Call  routine  to  calculate 
mesh  time  cycle(s) 


Boundary  condition 
iteration  loop,  repeat 
until  convergence 


Loop  for  100 
time  steps 


Call  routine  to 
calculate  tooth 
pair  loads  for 
either  1  or  2 
stages 


Call  routine  to 
solve  the  dynamic 
equations  for 
current  time  step 


(continued) 


Figure  3:  General  Program  Flowchart 
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(continued) 


(continued) 


i£ure  3:  General  Program  Flowchart  (continued) 
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(continued) 


gure  3 :  General  Program  Flowchart  (continued) 


TCI 


Stage  1 


C.O  0.025  0.05 

time,  sec. 

TC1/TC2  =  0.05/0.025  =  2.0 

Number  of  Boundary  Condition  Cycles 
for  Steady  State  =  2 

TCI  =  Time  Cycle  for  stage  1 
TC2  =  Time  Cycle  for  stage  2 


time,  sec. 


TC2/TC1  =  0 . 0264/ . 00562  =4.7 


10  x  TC2  -  0.264  sec. 
47  x  TCI  =  0.264  sec. 


Number  of  Boundary  Condition 
Cycles  for  Steady  State  -  47 


Figure  4:  Two  Stage  Time  Cycle/Boundary  Condition  Diagram 
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converged 


E:r  ;:^p 

■  value 

-i::.  I;:.;.!. 

V... 

.  .  i 

4  6  8 

Iteration  Number 

1 

converged  / 
;  value 


4  6  8 

Iteration  Number 


Figure  5:  Sun  Gear  Displacement  and  Velocity  Boundary 
Conditions  vs.  Number  of  Iterations 
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Figure  7:  Displacement  Boundary  Conditions  vs. 

Number  of  Iterations 
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Figure  8:  Velocity  Boundary  Conditions  vs 
Number  of  Iterations 


V/i 


SINGLE 

STAGE 

CHECK  CASE,  OH-58 

1 

1  . 

1  . 

o 

8.8571 

24.6 

3 

b  . 

2  7 

35. 

99. 

5 

9  p 

1 . 375 

1.2 

1. 

3.0 

1. 

4 

1<*. 

5000  . 

500. 

500  . 

1 .  0 

5 

28  . 

.  026 

.2000 

.00 

.  0160 

.016 

5 

33. 

.  016 

O 

150. 

.01 

20. 

1 

651  . 

1. 

0-1 . 


Figure  9a: 


Single  Stage  OH-58  Test  Case  Input 
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Figure  10:  Two  Stage  OH-58  Input  Data  Set 
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Figure  11a:  Dynamic  Loads  for  OH-58  First  Stage 

Sun-Planet  Mesh 
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Figure  lib:  Dynamic  Loads  for  OH-58  First  Stage 
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Figure  13:  Procedure  to  Locate  Source  of  Instability 
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